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1.  Introduction  and  Overview 

The  goal  of  our  research  is  to  be  able  to  monitor,  understand,  and  ultimately  forecast 
the  behavior  of  the  space  plasma  regimes  surroimding  the  Earth.  The  first  step  in  the 
procedure,  accomplished  in  1988-1991,  was  a  series  of  papers  on  the  identification  of  the 
sources  (such  as  the  shocked  solar  wind)  of  the  dayside  regions  of  precipitation  into  the 
ionosphere  as  accurately  as  possible  from  the  DMSP  series  low-altitude  satellites.  The 
identification  were  then  automated  [Newell  et  ah,  1991],  allowing  massive  statistical 
studies.  The  second  step,  funded  by  this  AFOSR  grant  and  carried  out  in  1992-1995 
using  the  capabilities  created  by  the  first  step,  has  resulted  in  the  first  global  mapping  of 
these  space  plasma  regions  --  the  magnetosphere  --  into  the  ionosphere  of  the  Earth, 
under  a  variety  of  conditions.  For  example,  Newell  and  Meng  [1994]  showed  the 
dramatic  effects  solar  wind  pressure  has  on  the  mapping  of  plasma  regions  into  the 
ionosphere. 

The  third  step  continuing  in  the  same  logical  progression  is  the  identification  and 
monitoring  of  the  magnetosphere  through  multiwavelength  UV  auroral  imagery. 
However  this  requires  calibrating  the  UV  signatures  to  the  particle  specifications  of  the 
plasma  source  regions.  It  is  fairly  clear  that  the  UV  images  contain  the  information 
needed  for  global  monitoring  of  the  magnetosphere  in  principle,  but  it  is  first  necessary 
to  develop  a  methodology  for  deconvolving  the  plasma  distribution  and  behavior  from 
the  images.  Some  very  preliminary  efforts  using  Polar  BEAR  and  DMSP  satellites  have 
already  conducted  [Newell  et  ah,  1992],  but  unfortunately  Polar  BEAR  had  only  limited 
capabilities  in  terms  of  choosing  the  wavelengths  that  could  be  simultaneously 
monitored,  and  it  is  not  possible  to  accurately  determine  the  magnetosphere's  plasma 
distribution  with  this  satellite. 

The  future  DMSP  satellites  will  contain  the  SSUSI  imager  which  will  allow  a  more 
fruitful  combination  of  wavelengths.  In  anticipation,  Meng  [1994]  has  carried  out  some 
of  the  preliminary  analysis  necessary  to  exploit  the  new  data  set,  by  converting  the 
particle  maps  of  plasma  distributions  [Newell  and  Meng,  1992]  into  intensities  at 
various  wavelengths. 

In  addition  to  the  logical  chain  above,  leading  to  the  monitoring  of  the 
magnetosphere  from  global  auroral  imagery,  a  second  promising  path  has  emerged 
from  work  on  the  nightside  precipitation.  It  turns  out  that  in  many  ways  the  state  of 
the  magnetosphere  can  be  more  accurately  determined  from  a  careful  study  of  nightside 
precipitation  than  anyone  had  supposed.  The  results  of  our  research  into  this  problem 
over  the  last  1.5  years  are  quite  significant,  and  we  now  devote  some  space  to  explaining 
them. 

Our  investigation  into  systematizing  nightside  precipitation  centered  on  extracting 
the  features  of  greatest  geophysical  significance.  The  proposed  system,  which  includes 
operational  definitions  and  has  been  automated,  consists  of:  boundary  (bl)  The  "zero- 
energy"  convection  boundary  (often  the  plasmapause);  boundary  (2e)  The  point  where 

the  large-scale  gradient  of  electron  average  energy  with  latitude,  dEe/dA.,  switches  from 
positive  to  <  0  (the  start  of  the  main  plasma  sheet);  boimdary  (2i)  The  ion  high-energy 
precipitation  cutoff  (the  ion  isotropy  boimdary  or  the  start  of  the  tail  current  sheet); 
boundary  (3a,b)  The  most  equatorward  and  poleward  electron  acceleration  events 
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(spectra  with  "monoenergetic  peaks")  above  0.25  ergs/cm^  s;  boundary  (4s)  Transition 
of  electron  precipitation  from  unstructured  on  a  >10  km  spatial  scale  (spectra  have  0.6- 
0.95  correlation  coefficients  with  neighbors)  to  structured  precipitation  (correlation 
coefficient  usually  0.4  and  below);  boimdary  (5)  The  poleward  edge  of  the  main  auroral 
oval,  marked  by  a  spatially  sharp  dropoff  in  energy  fluxes  by  a  factor  of  at  least  4  to 
levels  below  that  typical  of  the  auroral  oval;  boundary  (6)  The  poleward  edge  of  the 
subvisual  drizzle  often  observed  poleward  of  the  auroral  oval. 

As  has  been  mentioned,  the  boundaries  were  chosen  so  that  they  (1)  contained  as 
much  geophysical  significance  as  possible,  and  (2)  were  amenable  to  automated 
identification.  What  follows  in  section  2  is  a  more  detailed  explication  of  each 
boundary  with  respect  to  geophysical  significance;  in  section  3  a  brief  summary  of 
operational  identificatioris  will  be  give. 

2.  The  Structure  of  Nightside  Precipitation:  Physical  Significance  of  the  Boundaries 

As  in  many  situations,  an  example  is  clearer  than  words  alone;  hence  consider  Figure 
1,  which  presents  a  DMSP  F7  pass  at  11:23  UT  on  January  3,  1985.  We  consider  in  turn 
each  of  the  most  geophysically  significant  boundaries  that  can  be  identified  in  this 
picture.  Because  operational  definitions  which  are  robust  require  an  attention  to  detail 
that  is  tedious  to  many,  the  present  section  gives  conceptually  oriented  definitions, 
with  the  fine  print  reserved  for  Section  3. 

(1)  The  "zero-energy"  convection  boimdary  (11:27:13  UT).  Zero-energy  electrons  and 
ions  have  no  curvature  and  gradient  drifts,  hence  they  should  share  a  common 
equatorward  boundary;  one  which  is  determined  purely  by  the  electric  and  magnetic 
field  configuration.  (That  is,  the  boundary  which  results  only  from  a  consideration  of 
ExB  drift  effects).  The  DMSP  low-energy  ion  detector  has  an  extremely  large  geometric 
factor,  which  makes  it  possible  to  observe  such  coincidences  despite  the  comparatively 
low  fluxes  of  ions  at  low  energies.  Newell  and  Meng  [1988b]  have  reported  that  in  the 
dusk  and  midnight  sectors  the  electron  and  ion  zero  energy  cutoffs  indeed  coincide  on 
80%  of  the  passes.  To  maintain  operational  unambiguity,  we  propose  that  the  zero- 
energy  electron  and  ion  boundaries  be  separately  defined  (denoted  ble  and  bli).  Then, 

when  the  two  boundaries  indeed  coincide  to  within  0.25°  magnetic  latitude,  one  may 
reasonably  say  that  a  zero-energy  convection  boundary  exists.  Any  model  electric  and 
magnetic  field  of  the  magnetotail  specifies  a  unique  position  for  this  boundary  (as  a 
function  of  MLT);  hence  observation  of  the  boundary  provides  a  direct  comparison 
between  theory  and  reality.  In  some  theoretical  formulations,  the  zero-energy 
convection  boundary  is  also  the  plasmapause  location  [Nishida,  1966],  and  indeed 
observations  of  electron  data  support  this  association  [Galperin  et  al.,  1977;  Horwitz  et 
al.,  1986;  Sauvaud  et  al.,  1983].  However  the  zero-energy  boundary  observed  at  any 
given  time  does  not  necessarily  represent  a  steady-state  convection  boundary,  since 
magnetotail  convection  is  highly  dynamic  [Mauk  and  Meng,  1983].  Indeed,  the  low- 
energy  equatorward  ion  precipitation  often  contains  plasma  originally  of  ionospheric 
origin,  apparently  injected  at  high  latitudes  often  in  association  with  auroral  arcs 
[Bosqued  et  al.,  1986].  However  its  subsequent  convection  towards  lower  latitudes  is  the 
result  of  ExB  drift  effects.  To  maintain  operational  unambiguity,  and  because  the  zero- 
energy  boundary  has  theoretical  importance  in  its  own  right,  we  do  not  here  stress  the 
connection  with  the  plasmapause. 
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(2e)  Hie  poleward  edge  of  the  dEe/d^  >  0  region  (11:25:47  UT).  It  has  long  been  known 
that  low-energy  electrons  in  the  plasma  sheet  reach  closer  to  Earth  than  do  higher 
energy  electrons  [Vasyliunas,  1968;  Schield  and  Frank,  1970;  Fairfield  and  Vinas,  1984]. 
The  higher  the  energy  of  the  electrons  measured,  the  further  from  the  Earth  they 
appear  to  have  a  cutoff  (some  exceptions  exist,  such  as  when  a  dispersionless  injection 
occurs).  As  a  low-altitude  spacecraft  moves  poleward  from  boundary  le,  progressively 

higher-energy  electrons  are  observed;  so  dEe/dX.  >  0  (where  Ee  is  the  e-  average  energy). 
As  one  reaches  the  main  plasma  sheet,  electrons  of  all  energies  are  observed.  Further 

poleward  the  overall  trend  is  for  dEe/ dX  <  0  (the  region  of  negative  gradient  has  a  slope 
of  smaller  magnitude  and  exhibits  more  fluctuations  than  does  the  region  of  positive 
gradient).  This  is  simply  because  the  plasma  sheet  is  progressively  colder  further  from 

the  Earth.  The  point  where  dEe/ dX.  =  0  is  one  measure  of  the  start  of  the  main  plasma 
sheet  -  or  in  the  terminology  of  Feldstein  and  Galperin  [1985]  (hereinafter  FG85)  the 
true  start  of  the  central  plasma  sheet. 

(2i)  The  high-energy  ion  equatorward  precipitation  cutoff  (precipitating  energy  flux 
maximum  -  11:26:00  UT).  This  boundary  is  also  the  Isotropy  Boundary  (IB)  of  Sergeev 
et  al.  [1983].  It  is  probably  the  best  and  most  direct  proxy  for  the  location  of  the 
Earthward  edge  of  the  current  sheet.  Consider  ions  in  the  energy  range  from  a  few  keV 
to  tens  of  keV  (30  keV  for  DMSP).  Ions  in  this  energy  range  increase  in  temperature 
and  energy  flux  with  declining  latitude,  apparently  as  a  result  of  adiabatic  acceleration 
as  plasma  convects  Earthward  in  the  magnetotail  [Galperin  et  al.,  1977].  This  steady 
temperature  increase  terminates  with  a  relatively  sharp  equatorward  precipitation 
cutoff.  However  in  the  high-altitude  inner  magnetosphere,  ions  do  not  disappear  at 
the  L-shell  value  of  the  high-energy  ion  precipitation  boimdary  [e.g.,  Lui  et  al.,  1987]. 
Instead  what  changes  is  that  the  ions  become  trapped,  and  cease  to  precipitate  in 
measurable  quantities.  Poleward  of  the  precipitation  boimdary  at  any  particular  energy, 
the  ions  are  highly  isotropic  [Bernstein  et  al.,  1974;  Sergeev  et  al.,  1993].  It  has  thus  been 
suggested,  and  even  sucessfully  modeled  in  some  detail,  that  the  ions  maintain  their 
isotropy  by  pitch  angle  scattering  in  the  tail  current  sheet  [Lyons  and  Speiser,].  The 
physical  mechanism  is  quite  simple:  ions  cannot  maintain  pitch  angle  while  bending 
around  field  lines  which  have  a  radius  of  curvature  comparable  to  the  ion  gyroradi 
[Sergeev  et  al.,  1983].  This  explanation  also  accounts  for  the  dispersion  in  the  high- 
energy  ion  cutoffs  [Sergeev  et  al.,  1993].  The  larger  gyroradii  of  higher  energy  ions 
means  that  they  scatter  off  field  lines  with  smaller  radii  of  curvature  than  do  the  lower 
energy  ions;  hence  the  higher-energy  ions  maintain  isotropy  further  Earthward. 

Neither  the  tail  current  sheet  nor  the  precipitating  high-energy  ions  have  a  sharply 
defined  boundary.  Operationally  we  propose  to  use  the  ion  precipitating  energy  flux 
peak  (integrated  over  the  range  3  keV  to  30  keV),  which  universally  occurs  close  to  the 
equatorward  boimdary  of  the  high-energy  ion  precipitation,  as  the  location  of  b2i.  The 
geophysical  significance  of  the  boundary  is  that  it  represents  a  good  approximation  to 
the  Earthward  edge  of  the  tail  current  sheet.  Sergeev  and  Gvozdevsky  [1995]  have 
demonstrated  that  the  latitude  of  this  ion  isotropy  boundary  has  a  very  high  correlation 
(r~0.9)  with  the  magnetic  field  inclination  (degree  of  stretching)  measured 
simultaneously  at  the  geomagnetic  equator. 
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A  boundary  closely  related  conceptually  is  the  >  30-40  keV  electron  trapping  boundary 
for  electrons,  which  we  term  b2t.  Although  we  cannot  directly  identify  it  in  our  DMSP 
data  base,  it  has  a  long  history  of  being  considered  useful.  Thought  in  the  1960s  to 
represent  the  open/closed  field  line  boundary,  b2t  is  now  generally  recognized  as 
another  measure  of  where  field  lines  begin  to  be  significantly  stretched.  Because  of  the 
smaller  g5nroradii  of  electrons,  b2t  usually  lies  a  short  distance  poleward  of  b2i.  There  is 
some  reason  to  believe  that  the  stable  trapping  boimdary  approximately  corresponds  to 
the  Winningham  et  al.  [1975]  (hereinafter  WYAH75)  boimdary  between  CPS  and  BPS; 
and  indeed  Weiss  et  al.  [1992]  define  this  latter  boundary  by  the  trapping  boundary. 

(3a,b)  The  most  equatorward  and  poleward  e-  acceleration  events.  In  the  literature 
many  proxies  for  identifying  the  region  of  discrete  auroras  exist.  For  example  it  appears 
that  most  electron  acceleration  events,  and  certainly  those  of  high  accelerating  potential 
values,  occur  on  the  stretched  field  lines  that  lie  poleward  of  the  >  40  keV  e-  stable 
trapping  boundary  [Frank  and  Ackerson,  1971].  It  is  quite  feasible  to  individually 
examine  each  electron  spectrum  and  determine  whether  it  shows  signs  of  a  field- 
aligned  accelerating  potential.  Therefore  we  include  boundaries  3a  and  3b,  which,  based 
on  the  examination  of  each  individual  spectrum,  are  the  furthest  equatorward  and 
furthest  poleward  sites  of  electron  acceleration.  A  spectrum  is  identified  as  accelerated 
if  it  either  has  a  monoenergetic  peak,  or  a  sharp  cutoff  above  the  spectral  peak  (more 
detail  is  available  in  Newell  et  al.  [1996]).  Although  we  have  opinions  about  the  likely 
location  of  these  boundaries,  it  is  best  that  they  be  identified  separately  from  all  other 
precipitation  boundaries. 

Boundary  4s—  the  onset  of  spatial  structure  in  electron  precipitation  (on  a  scale  >  5-10 
km)  (11:25:14).  The  BPS/CPS  distinction,  no  matter  what  misconceptions  became 
associated  with  it,  has  persisted  primarily  because  many  nightside  crossings  do  seem  to 
have  a  highly  structured  region  and  a  relatively  unstructured  region  of  auroral  electron 
precipitation.  If  this  structured /unstructured  distinction  represents  something 
fundamental,  there  should  be  a  quantitative  way  of  making  the  distinction  within  the 
precipitation  data  itself,  i.e.,  one  which  does  not  depend  on  such  additional  factors  such 
as  the  boundary  of  the  radiation  belts  (b2t). 

To  move  from  a  qualitative  description  ("structured")  to  a  quantitative  description, 
we  investigated  the  behavior  of  the  correlation  coefficient  of  individual  spectra  with 
their  neighbors.  Figure  2  shows  a  plot  of  the  running  average  of  the  correlation 
coefficient  of  each  spectra  in  Figure  2  with  the  previous  5  spectra.  (Our  algorithm 
suppresses  spectra  with  fluxes  below  0.25  ergs/cm2  s  by  halving  the  correlation 
coefficient.)  Although  even  within  the  "BPS"  region  the  majority  of  the  individual 
spectra  do  not  show  evidence  of  field-aligned  acceleration.  Figure  2  shows  that  the 
entire  region  is  indeed  structured  in  the  sense  that  each  point  correlates  only  poorly 
with  its  neighbors.  We  therefore  introduce  boundary  4s,  the  structured  precipitation 
boundary,  defined  by  the  point  where  the  correlation  coefficient  drops  from  the  0.95- 
0.60  range  to  below  0.60.  Thus  b4s  marks  a  fundamental  change  in  the  character  of  the 
electron  precipitation,  from  a  spatially  unstructured  region  to  a  highly  structured 
region. 

Boundary  5(e,i)  —  The  poleward  boundary  of  the  main  auroral  oval  (11:24:36  UT).  The 
precipitating  energy  flux  in  the  auroral  oval  t5q)ically  drops  by  about  an  order  of 
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magnitude  over  a  short  distance  (usually  <~0.2o).  This  dramatic  dropoff  occurs  in  both 
the  electrons  and  ions,  although  not  always  in  precisely  the  same  location  for  the  two 
species.  In  active  times,  poleward  of  this  sharp  dropoff  there  is  often  only  a  narrow 
region  of  high  energy  (~10  keV)  but  diffuse  electrons  at  a  very  low  flux  levels.  (This 
faint  "overhang"  of  high-energy  electrons  seems  to  be  both  an  unremarked  feature  in 
the  literature  and  a  common  feature  in  nature).  Sometimes,  especially  in  quiet  times, 
including  quieting  times  following  a  substorm,  a  subvisual  drizzle  can  extend  poleward 
of  the  oval.  The  conceptual  definition  of  the  poleward  edge  of  the  main  auroral  oval  is 
that  the  precipitating  fluxes  drop  by  a  factor  of  at  least  4  over  a  short  distance  to  values 
below  3x1010  ev/cm2  s  str  (e-)  or  1010  ev/cm2  s  str  (ions).  We  emphasize  that  even  for 
northward  IMF  conditions  such  sharp  dropoffs  still  occur,  and  separate  the  main 
auroral  oval  from  the  polar  cap  precipitation  (which,  to  be  sure,  usually  resembles  the 
structured  plasma  sheet  precipitation  although  not  at  the  intensity  of  the  main  oval). 
Boundary  6.  The  poleward  boundary  of  the  subvisual  drizzle  (11:23:53  UT).  Poleward 
of  the  main  auroral  oval  is  a  weak  subvisual  drizzle  that  differs  from  polar  rain  in 
several  ways.  The  subvisual  drizzle  usually  includes  weak  ion  as  well  as  electron 
precipitation;  the  drizzle  is  structured,  while  polar  rain  is  comparatively  homogeneous; 
and  the  typical  electron  energies  are  a  bit  higher  than  normal  for  polar  rain.  Polar  rain 
is  intense  poleward  of  the  dayside  oval,  and  declines  gradually  in  intensity  moving 
towards  the  nightside  oval.  The  subvisual  drizzle  extends  poleward  from  the  nightside 
oval  and  terminates  either  when  fluxes  drop  to  background  levels,  or  (less  commonly) 
when  a  smooth  coherent  electron  polar  rain  signature  is  encoimtered. 

3.  Operational  Definitions  of  Nightside  Precipitation 

In  this  section  we  give  the  (sometimes  tedious)  detail  required  to  operationally  define 
the  various  boundaries  unambiguously.  The  integral  parameters  are  denoted  n,  jE,  and 
E,  referring  to  density,  energy  flux,  and  average  energy  respectively.  These  quantities 
can  be  subscripted  by  "e"  and  "i"  for  electrons  and  ions.  Better  results  are  obtained  if  the 
terms  jep(El,E2)  etc  are  introduced,  referring  to  the  "partial"  electron  energy  flux 
between  El  and  E2.  Strictly  speaking  all  instruments  measure  only  a  partial  values  of 
"integral"  parameters  between  the  upper  and  lower  ranges  of  their  detector  (some 
confusion  in  the  literature  exists  from  the  neglect  of  this  point,  as  in  the  case  of  the  cusp 
electron  average  energy).  As  a  practical  consideration  here,  such  partials  improve 
greatly  the  identification  of  boundaries  while  also  making  the  results  somewhat  less 
dependent  on  the  specific  detector  (SSJ4)  that  they  were  designed  for.  The  units  used 
below  are  eV  for  E,  and  loglO  eV/cm2  s  str  for  jep,  jip. 

For  each  boundary,  we  start  by  giving  the  rule  which  works  most  of  the  time,  and 
follow  with  caveats  and  more  detail  on  handling  special  cases. 

—  ble,i  (zero-energy).  The  algorithm  moves  from  lower  latitudes  to  higher,  comparing 
the  average  of  jip(El,E2)  and  jep(El,E2)  (ordinarily  El  and  E2  are  the  2  lowest  channels) 
over  the  3  previous  spectra  with  the  3  succeeding  spectra.  An  increase  in  jip  by  a  factor 
of  2  marks  the  onset  of  the  zero  energy  boundary,  which  is  separately  determined  for 
the  two  species.  This  jump  is  significant  only  if  it  also  significantly  exceeds  the 
background  counts  gotten  by  averaging  over  several  equatorward  seconds.  If  jep  rises  to 
a  value  above  8.0  (jip  reaches  6.5)  a  factor  of  only  1.6  jump  is  acceptable  in  determining 
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the  zero-energy  boundary.  If  jep  >  8.25  (jip>6.9),  it  is  assumed  that  the  boundary  has 
been  reached,  even  if  no  jump  in  the  value  of  the  fluxes  is  measured. 

Special  cases:  The  energy  range  considered  (El  to  E2)  in  the  partials  depends  on 
whether  photoelectrons  are  present,  and  whether  the  spacecraft  is  charged  to  -28  V.  The 
former  can  be  identified  by  a  sharp  dropoff  in  electron  fluxes  above  68  eV  at  latitudes 
below  the  auroral  zone;  the  latter  by  a  sharp  cutoff  above  the  32  eV  ion  channel.  In  the 
absence  of  these  effects,  the  channels  are  set  to  the  lowest  available,  El=32,  E2=47.  If  the 
spacecraft  is  charged  to  -28  V,  the  ion  channels  are  set  to  El=47,  E2=68.  If  photoelectrons 
are  present  (rare  on  the  nightside),  the  next  available  "clean"  channels  are  100  and  145 
eV.  Finally  sometimes  isolated  noise  can  cause  false  positives,  as  by  radiation  belts 
counts.  Thus  a  "checkble"  routine  performs  a  double  check  by  simply  examining  the 
next  several  seconds.  If  instead  of  rising  fluxes  as  the  auroral  oval  is  entered,  the  next 
few  seconds  exhibit  a  dropoff  in  fluxes,  the  identification  of  ble  is  inaccurate  and  the 
search  resumes  towards  increasing  latitudes. 

—  b2e  (plasma  sheet  start).  This  algorithm  locates  the  first  point  poleward  of  ble  where 
dEe/dl  <0.  A  sliding  3-s  average  Ee  is  compared  to  that  of  any  three  consecutive 
seconds  over  the  next  9  s.  If  the  value  of  Ee  does  not  rise  over  this  interval,  the 
boundary  b2e  has  been  located. 

Special  cases:  If  the  correct  boundary  b2e  has  been  found,  one  expects  to  have  entered 
the  main  plasma  sheet.  Hence  if  jEe  <  11  (or  if  both  jEe  <  11.5  and  Ee  <1000  eV)  then 
the  next  30  seconds  are  examined,  to  see  if  a  spectra  exists  with  both  higher  Ee  and  a  je 
larger  by  at  least  0.3  (i.e.,  a  factor  of  2  difference  in  the  energy  flux).  This  double 
checking  is  done  only  under  these  suspicious  circumstances  of  low  fluxes,  because  it 
risks  encountering  a  point  with  high  Ee  and  high  je  simply  because  strong  field-aligned 
acceleration  is  encountered. 

~b2i  (ion  isotropy  boxmdary).  This  boundary  is  defined  by  the  precipitation  flux 
maxima  for  ions  3  keV  and  above.  The  high-energy  ions  behave  less  variably  than  do 
the  electrons,  so  fewer  precautions  are  needed.  Thus  a  sliding  average  of  jip(3  keV,30 
keV)  over  2  seconds  is  compared  to  jip  for  any  three  consecutive  seconds  over  the  next 
10  s  further  poleward  to  determine  whether  the  maxima  has  been  located.  The  smallest 
maxima  acceptable  is  10.5. 

Special  cases:  Sometimes  "nose"  events  occur,  namely  the  injection  of  high  energy  ion 
regions  slightly  detached  from  the  rest  of  the  auroral  oval  [Konradi  et  al.,  1975;  Sanchez 
et  al.,  1993].  Such  events  are  identified  as  a  local  maxima  followed  by  a  decline  and 
subsequent  rise  to  a  global  energy  flux  maxima.  The  detached  (or  partially  detached) 
maxima  is  discarded. 

Additional  Note:  It  is  possible  (in  times  of  prolonged  quiet)  for  b2i  to  occur 
equatorward  of  bli;  but  by  definition  b2e  must  always  lie  poleward  of  ble. 

— ^b3a,b  (most  equatorward  and  poleward  e-  acceleration  events).  Each  individual 
spectra  is  examined  for  evidence  of  a  monoenergetic  peak.  This  can  either  be  a  single 
channel  with  a  differential  energy  flux  5  times  larger  than  any  other;  or  a  sharp  drop  by 
at  least  a  factor  of  10  below  the  electron  differential  energy  flux  peak.  Details  on  this 
algorithm,  including  special  cases,  are  presented  in  a  separate  paper  [Newell  et  al.,  1996]. 
The  most  equatorward  and  poleward  individual  spectra  showing  such 
"monoenergetic"  peaks  are  flagged  as  boundaries  b3a  and  b3b,  without  regard  to  the 
other  boundary  locations. 


7 


— ^b4s  (structured /unstructured  boundary).  The  counts  in  the  various  channels  for  a 
given  spectrum  are  correlated  with  the  corresponding  counts  in  the  5  previous  spectra, 
and  the  5  resulting  correlation  coefficients  are  averaged  (<r>).  By  definition  b4s  lies 
poleward  of  b2e  and  b2i.  When  the  sum  of  seven  consecutive  <r>  drops  below  4.0,  the 
search  is  halted.  b4s  is  set  to  be  the  furthest  poleward  spectra  within  the  current  group 
of  7  seconds  which  has  <r>  >  0.60. 

Special  cases:  If  the  energy  fluxes  result  only  in  an  aurora  which  is  subvisual  (je<10.7, 
or  0.25  ergs/ cm2  s)  the  correlation  coefficient  is  suppressed  (halved);  hence  low  flux  but 
homogeneous  regions  such  as  polar  rain  are  automatically  excluded. 

— ^b5e,b5i  (poleward  edge  of  main  oval).  The  boundaries  are  computed  separately,  but 
using  the  same  procedure.  An  average  of  jE  for  the  previous  12  seconds  is  compared 
with  jE  for  the  succeeding  12  seconds.  When  a  dropoff  of  a  factor  of  4  is  located,  a 
provisional  b5  boundary  is  determined.  Note  that  this  algorithm  emphasizes  locating  a 
sharp  drop  in  the  flux  levels. 

Special  cases.  The  next  30  s  are  double  checked  (35  s  for  e-)  to  make  sure  the  dropoff 
remained  below  auroral  energy  fluxes.  If  the  (log)  average  jE  has  not  dropped  below 
about  9.7  for  ions  or  10.5  for  electrons,  the  search  continues  for  the  corresponding  b5 
boundary.  Such  a  large  search  ahead  is  needed  because  of  features  such  as  the  "double 
oval"  [Elphinstone  et  al.,  1995]  in  which  fluxes  temporarily  drop  below  oval  levels  only 
to  rise  again  clearly  to  oval  values. 

— ^b6.  The  poleward  boundary  of  the  subvisual  drizzle  is  defined  by  the  point  where 
either  polar  rain  is  encountered  (identified  by  the  presence  of  unstructured  electrons 
and  no  ions)  or  je  drops  below  10.4  and  ji  drops  below  9.6. 

Special  cases:  Tlie  drizzle  is  defined  by  weak  structured  fluxes  with  ions  and  electrons 
above  noise  levels.  By  checking  the  computed  average  energy  one  can  infer  whether 
noise  is  significant:  For  example  if  coimts  are  randomly  distributed  across  all  channels 
one  would  obtain  E  =  15  keV  for  DMSP.  If  Ee  <  500  eV,  even  a  flux  as  low  as  10.0  is 
acceptable.  The  lower  the  average  energy  the  lower  a  flux  level  is  acceptable  as  not 
representing  noise. 

4.  Summary 

The  major  progress  made  within  the  4  year  cycle  (including  the  year  under  a  no-cost 
extension)  of  AFOSR  grant  F49620-92-J-0196  occurred  in  two  major  areas.  The  first  is  a 
further  development  of  the  research  line  originating  in  our  multi-year  effort  to  identify 
the  source  region  of  dayside  precipitation  regions  as  accurately  as  possible.  Because  the 
ionosphere  is  a  much  more  compact  and  cheaply  accessible  region  to  sample  than  is  the 
magnetosphere  to  which  it  is  geomagnetically  connected,  the  benefits  of  this  approach 
to  investigating  and  monitoring  the  state  of  the  magnetosphere  are  manifold. 

Major  research  accomplishments  along  this  line  of  endeavor  include  creating  maps  of 
the  ionosphere's  connection  to  the  magnetosphere  under  varying  conditions  of  solar 
wind  pressure  [Newell  and  Meng,  1994].  Since  the  logical  extension  of  this  work  for 
monitoring  the  magnetosphere  includes  the  identification  of  magnetospheric 
boundaries  from  multi-spectral  imagery,  some  preliminary  feasibility  work  in  this 
direction  has  been  carried  out  [Newell  et  al.,  1992;  Meng,  1994]. 

Because  much  of  magnetospheric  d5mamics,  and  in  particular  large  magnetic  storms 
affecting  Air  Force  operational  environments,  occurs  on  the  nightside,  we  have  within 
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the  last  grant  cycle  considered  the  problem  of  monitoring  of  the  state  of  the 
magnetosphere  from  nightside  precipitation  structure.  The  boundaries  of  greatest 
geophysical  significance  were  identified,  and  automated  algorithms  created  and  tested 
for  their  identification.  Although  our  system  of  dayside  identifications  has  had 
considerable  academic  impact,  being  both  widely  cited  and  widely  incorporated  by  other 
researchers  into  their  investigations  of  the  magnetosphere  from  groimd  based  and  low- 
altitude  measurements,  it  has  not  directly  led  to  any  immediate  improvements  in  space 
weather  forecasting.  Our  new  work,  completed  in  the  latest  grant  cycle  and  discussed  in 
sections  2  and  3  of  the  present  report,  extends  the  sophisticated  monitoring  of  the 
magnetosphere  into  the  dynamic  region  of  the  nightside.  We  expect  that  the  impact  on 
other  researchers  of  the  near-Earth  space  system  will  be  comparable  to  that  of  the 
dayside  work;  moreover  we  believe  that  the  new  system  will  be  more  likely  to  direct 
lead  to  improvements  in  space  weather  forecasting. 
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